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A stochastic version of an inviscid dyadic model of turbulence, 
with multiplicative noise, is proved to exhibit energy dissipation in 
spite of the formal energy conservation. As a consequence, global 
regular solutions cannot exist. After some reductions, the main tool 
Q^ \ is the escape bahavior at infinity of a certain birth and death pro- 

Pi, ■ cess. 

^' 

+-i ■ 1. Introduction. The dyadic model of turbulence has been introduced in, 

among others, [15, 16, 27] and [13], as a simplified model of fluid dynamics 

equations in order to investigate a number of properties which are out of 

reach at present for more realistic models. In this paper we study a suitable 

\^ , random perturbation of the classical dyadic model under which we are able 

tH- I to prove anomalous dissipation of energy. 

O ■ On a complete filtered probability space {Q,Ft,P), let {Wn)n>i be a se- 

^^ I quence of independent Brownian motions. Consider the infinite system of 

stochastic differential equations in Stratonovich form 

O . (1) dXn = (A;„„iX„„i - knXnXn+l) dt + kn^iXn^l o dWn-l - KXn+l o dWn 



for n > 1, with Xo(i) = 0. Denote by P the Hilbert space of real square 
summable sequences x = {xn)n>i and set ||3;||^ = Yl'^=i^'n- We call energy 
of X{t) := {Xn{t))n>i the quantity £{t) := ^||X(t)||^. Assume for simplicity 
to have a deterministic initial condition 

X{0)=X, X={Xn)n>l^l'^- 

The sequence of positive numbers (A:n)n>i will be specified later on; the most 
natural case in analogy with fluid dynamics is kn = A" for some A > 1 . 
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2 D. BARBATO, F. FLANDOLI AND F. MORANDIN 

System (1) is formally energy preserving. By the Stratonovich form of 
Ito's formula (see [19]), we have 

ldXl = XnOdXn 

= {kn~lX.^_iXn — knX^Xn+l) dt 

+ kn-lXn-lXn o dWn-1 - A;„,X„X„+i O dWn- 

If we sum formally these identities and use the boundary condition Xo(t) = 0, 
we readily have ^d^^^i X^ = 0, namely, 

£■(*)= £"(0), P-a.s. 

The aim of this paper is to prove rigorously an opposite statement, a prop- 
erty that we could call anomalous dissipation. We need the notion of energy 
controlled solution that will be given in the next section. 

Theorem 1. Assume kn = A" for some A > 1. Given x G /^, let X{t) be 
the unique energy controlled solution of equation (1). Then, for all t > 0, 

p{£{t)=£m<i 

and for all e > there exist t such that 

P{£{t) < e) > 0. 

Moreover, if £{0) is sufficiently small, then £{t) decays to zero at least ex- 
ponentially fast both almost surely and in C} . 

As a consequence of this theorem, in Section 7 we will also prove that 
global regular solutions cannot exist. 

The proof of the theorem is built along Sections 3-6 and concluded in 
Section 6.1. 

Being inspired by fluid dynamics, the results of Theorem 1 could be in- 
terpreted as a form of turbulent dissipation. Dynamically speaking, it is 
a dissipation due to a very fast cascade mechanism; energy moves faster and 
faster from low to high wave numbers n and escapes to infinity in finite time. 

Results of anomalous dissipation for linear stochastic systems with addi- 
tive noise have been proved by [22, 23]. The notion of anomalous dissipation 
of these papers is different and based on invariant measures, but conceptu- 
ally the question is the same. The dyadic model of the present paper is, to 
our knowledge, the first nonlinear stochastic case where anomalous dissipa- 
tion is proved to occur. Moreover, the proof is entirely different from those 
of the additive noise linear stochastic case. 

Nonlinear models with anomalous dissipation have been discovered before 
in the deterministic case (see [3, 5, 7-9, 15-17, 27]). Our model is a mul- 
tiplicative random perturbation of models of these forms. However, let us 
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stress that the proof given here is totahy different from the proofs of the de- 
terministic hterature, based on monotonicity and positivity properties of the 
deterministic part of system (1). These properties are lost in the stochastic 
case. 

Theorem 1 remains true when A;„ < Cn" for some q > 1. However, we 
make here the assumption fc„ = A" in analogy with the deterministic litera- 
ture on dyadic models. 

The proof is based on three main ingredients: (i) Girsanov's transform 
allows us to reduce the problem to a linear stochastic equation; (ii) second 
moments of components satisfy a closed system, without moments of prod- 
ucts; (iii) this closed system is the forward equation of a birth and death 
process, the escape of which at infinity can be understood. 

Some issues in this procedure are not trivial. One of them is the equiv- 
alence of laws on infinite time horizon (see Proposition 18). Its proof is 
nonstandard; moreover, it is restricted to a range of values of parameters, 
the generalization being open. Concerning the idea that square moments 
could satisfy a closed equation, related to jump process, we have been in- 
spired by previous works ([1, 10, 11] and [21]), however, devoted to different 
models; the link here with the nonlinear model and the stochasticity is more 
transparent via Girsanov's transform. 

1.1. The multiplicative noise in Euler's equations. The noise we intro- 
duced in equation (1), motivated by energy conservation, may appear pecu- 
liar from the physical point of view. Nevertheless, it is the natural choice if 
we compare the dyadic model with some other equations of fluid dynamics 
like Euler's equations or diffusions of passive scalars. 

Let us give a picture of this analogy in the case of Euler's equations, which 
have the form 

du 

— — h u ■ Vu + \7p = 0, div n = 

ot 

with appropriate initial and boundary conditions {u and p are the veloc- 
ity and pressure field, resp.). Let us think of periodic boundary conditions 
for sake of simplicity. The Lagrangian motion of particles is given by the 
equation 

^ = «(.,y(*))^ 

A natural way to randomly perturb Euler dynamics (see [24]) is by adding 
a white noise to the Lagrangian motion, 

dY{t) = u{t,Y{t)) dt + ^aj{t,Y{t)) dW^ (t), 
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where {W^ {t))t>o are independent Brownian motions and ctj are given vector 
fields. By standard rules of stochastic calculus, applied formally, one can see 
that Euler equations take the stochastic form (with a new pressure p) 

(2) du+lu-Vu + Vpldt + '^aj-Vuo dW^{t) = 0, divn = 0, 

j 

where Stratonovich operation has to be used. Rigorous results and physi- 
cal arguments in support of this kind of stochastic perturbation of the La- 
grangian motion and the corresponding PDE with multiplicative Stratonovich 
noise (in the viscous case) can be found in [24, 25]. In addition, let us men- 
tion the wide literature on stochastic passive scalar equations (see, e.g., [20]) 
where multiplicative Stratonovich noise of the form above is used. 
In abstract form, equation (2) takes the form 

(3) du + B{u,u)dt + B{odW,u), 

where W{t) := X],crj(2;)VF-'(t) and B{u,v) := u ■ Vv. Notice that for suffi- 
cient regular u and v the following identity holds {B{u,v),v) = 0. 

In [8] the authors argued that after Fourier or wavelets transforms and 
certain simplifications, the deterministic system 

(4) — ^ = A;„_iX^,_i - A;„X„X„+i, n > 1, 

describes some idealized features of the deterministic equation -^ + B{u,u) = 
0. Equation (4) has the form -^X = B{X,X) where 



B{X,Y)(^n) — kn~lXn~lYn^i — knXnYn+l- 

Formally we have {B{X, Y),Y) = 0. Thanks to {B{u, u),u)=0 and {B{X, X), 
X) = the perturbation of u and X does not modify the energy balance (at 
a formal level). 

Standing this idealized discretization B{X, Y) of B{u, v), the natural ana- 
log of equation (3) is 

dXn + B{X,X)^^)dt + B{odW,X)^^), n>l, 

which is precisely system (1). 

2. Ito's formulation. For the rigorous formulation of equation (1) and 
a basic theorem of existence and uniqueness, we follow [4]. The Ito form of 
equation (1) is 

dXn = {kn-lXl^^ - knXnXn+l) dt + kn-lXn-1 dWn-1 
(5) 

- knXn+l dWn - \{kl + kl_^)Xn dt. 
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Let us define the concept of weak solution for tliis equation. By a filtered 
probability space {Q,Ft,P) we mean a probability space {Q,Fao,P) and 
a right-continuous filtration {Ft)t>o such that F^o is the u-algebra generated 

Definition 2. Given x € /^, a weak solution of equation (1) in P is 
a filtered probability space {Q,Ft,P), a sequence of independent Brownian 
motions {Wn)n>i on (fl,Fj,P) and an /^-valued stochastic process (X„)„>i 
on {Q,Ft,P) with continuous adapted components X„, such that 

Xnit) =Xn+ / (kn-^lXl^iis) - /c„X„(s)X„+i (s)) ds 

Jo 

+ / kn~iXn^l{s) dWn~l{s) - /c„X„+i(s) (iW„(s) 

JO JO 



j ^ikl + kl_^)Xnis)ds 



'0 

for each n > 1, with Xq = 0. We denote this solution by 

{n,Ft,p,w,x) 

or simply by X. 

Definition 3. We call energy controlled solutions the solutions of Def- 
inition 2 which satisfy 

(oo oo \ 

Y^Xl{t)<Y^xl]=l 
n=l n=l / 

for ah t > 0. 

The following simple proposition (proved in [4]) clarifies that a process 
satisfying (5) rigorously satisfies also (1). 

Proposition 4. If X is a weak solution, for every n > 1 the process 
{Xn{t))t>o is a continuous semimartingale, hence, the two Stratonovich in- 
tegrals 

ft f-t 

kn^iXn^i{s) o dWn-i{s) - / A;„X„+i(s) o dWnis) 
10 Jo 

are well defined and equal to 



t 1 rt 



kn-lXn-l{s) dWn-l{s) - - / k^_iXn{s) ds 



2./0 



- / knXn+l{s) dWnis) - - / /c^X„(s)ds. 

Jo ^ Jo 

Hence, X satisfies the Stratonovich equations (1). 
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The main result proved in [4] is the well posedness in the weak proba- 
bilistic sense in the class of energy controlled solutions. 

Theorem 5. Given (a;„) S P, there exists one and only one energy con- 
trolled solution of equation (1). 

3. Girsanov's transformation. Formally, let us write equation (5) in the 
form 

dXn = kn^iXn^iiXn-i dt + dWn~i) - A;„X„+i(X„ dt + dWn) 

~ 2\ n ~^ l^n—ll-^n dt. 

The simple idea is that Xn dt + dWn is a Brownian motion with respect to 
a new measure Q on (0, F), simultaneously for every n, hence, the equations 
become linear SDEs under Q. We use details about Girsanov's theorem 
that can be found in [26], Chapter VIII, and an infinite dimensional version 
proved in [6, 12, 18]. 

Assume that {Xn)n>i is an energy controlled solution. Due to the bound- 

edness of XlJ^i ^nif) [see (6)], the process It := - Yl'^=i /o ^n{s) dWn{s) is 
well defined, is a martingale and its quadratic variation [y, Y]t is 
Jp^^-^X^(s) ds. For the same reason, Novikov criterium applies, so 
N{Y)t := exp(yt — [y, y]t) is a strictly positive martingale. Define the set 
function Q on IJt>o -^* ^^ setting 

*l =A/-(F). = expUf;|'x„(s)<i»'„(s)-ljr'f;X,2(.)<!5J 



<^> dp 



n=l''^ •''J n=l 

for every t > 0. We also denote by Q its extension to the terminal ir-field F^. 
In general we cannot prove it is absolutely continuous with respect to P, 
but we shall see at least a case when this is true. Notice also that Q and P 
are equivalent on each Ft, by the strict positivity. Define 

Bn{t)=Wn{t)+ I Xn{s)ds. 

Jo 
Under Q, {Bn(t))n>i^te[o,T] is a sequence of independent Brownian motions. 
Since 

/ kn-lXn-i{s)dBn-l{s)= j kn-lXn-l{s) dWn-l{s) 

Jo Jo 

+ / kn~lXn^l{s)Xn-lis) ds 

Jo 
and similarly for J^ knXn+i{s)dBn{s), we see that 

Xn{t)=Xn{Q)+ I kn^iXn^l{s)dBn^l{s)- I /c„X„+i(s) dB„(s) 

Jo Jo 

/■* 1 

- / 2''kl + kl_^)Xn{s)ds. 
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This is a linear stochastic equation. Girsanov's transformation has removed 
the nonhnearity. Let us collect the previous facts. 

Theorem 6. // {Q,Ft,P,W,X) is an energy controlled solution of the 
nonlinear equation (1), then it satisfies the linear equation 

(8) dXn = K^iXn^i dBn-i - fc„X„+i dBn - |(A;^ + k'i-i)Xn dt, 
where the processes 



Bn{t) = Wn{t)+ [ Xn{s)ds 

Jo 



are a sequence of independent Brownian motions on {Q,Ft,Q), Q defined 
by (7). 

One may also check that 

dXn = kn-lXn-l o dBn-1 — knXn+1 o dBn 

so the previous computations could be described at the level of Stratonovich 
calculus. 

4. Closed equation for E^[Xl{t)]. Let {n,Ft,P,W,X) be an energy 
controlled solution of the nonlinear equation (1) with initial condition x gP 
and let Q be the measure given by Theorem 6. Denote by E*^ the mathe- 
matical expectation on {Q,Ft,Q). We have 

2 d-^n = -^n dXn + 2 d[Xn]t 

= ~2\kn + Ki-l)^ndi + dMn + 2(^71-1^^-1 + ^n^n+l)"*) 
2 '^■^n — ^-^n dXn + -Y-^n d[Xn\t 

= -\[kl + A;2_i)X^ dt + dMr, + 'iXl{kl_^Xl_, + klXl^,)dt, 
de^" = e^" dX^ + ie^' d{X^, X^], 
de^" = e^" dXn + ie^" d[Xn,Xn] 

= ■■■ e^"{kl + kl__^)Xndt + dMn + \e^^\{kl^iXl^^ + klxl^^) dt, 
where 

Mnit)= [ kn-lXn-lis)Xnis)dBn^lis) - [ A;„X„(s)X„+i(s) rf5„(s). 

Jo Jo 

Notice that 

(9) E^ / X^{t)dt< 00 

Jo 
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for each n > 1. Indeed, for an energy controlled solution, from (6) we have, 
with P-probability one, 

oo oo / oo \ '^ 

5]X^(i)<maxX2(t)J]x2(t)<K;4 . 

71=1 n=l \n=l / 

But P and Q are equivalent on Ft, hence, 

(oo / oo \ ^\ 

Y.<it)<\y.^i] 1=1- 

This implies (9). 

From (9), Mn{t) is a martingale for each n > 1. Moreover, 
'^^[Yl'^=i -^nit)] < C'o because {Xn)n>i is an energy controlled solution [again, 
as above, condition (6) is invariant under the change of measure P -H- Q 
on Ft] and thus, in particular, E*5[X^(t)] is finite for each n>l. From the 
previous equation we deduce the following. 

Proposition 7. For every energy controlled solution X, ¥fl[X'^{t)] is 
finite for each n > 1 and satisfies 

^E<^[Xl] = -{kl + kl_,)E<^[Xl] 

for t > 0. 

The first remarkable fact of this result is that E'5[X^] satisfies a closed 
equation. The second one is that this is the forward equation of a continuous- 
time Markov chain, as we shall discuss in the next section. See [1, 10] for 
different examples with the same structure. 

5. Associated birth and death process. In this section we will make thor- 
ough use of birth and death processes. We do not suppose that all the readers 
are familiar with the field, so we will be more detailed. 

Let us set 

Pnii) = Tr^^'^[X?i{t)], P{t) = {pn{t))n>l, t>0, 

and set also po{t) = 0. Introduce the positive numbers (Ar!,)n>i and (/^n)n>i) 
defined as 

\ — k"^ II — h^ 
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By Proposition 7, we have 
(10) < 



-lTPn(i) = -(An +/in)pn(i) + A„_ip„_i(t) + /i„+ip„+i(t), t > 0, 

Pn{0) ^' 



I 119 * 

We observe that X^^i J'n(i) = 1 when t = and, moreover, 

oo 

(11) ^Pnit)<l 

n=l 

for all t > (since X is an energy controlled solution). 

The system (10) can be conveniently put in matrix form, -^p{t) =p{t)A, 
where A is an infinite matrix with null row sums and nonnegative off- 
diagonal entries. 

In the theory of continuous-time Markov chains this is usually referred to 
as a qi-matrix. Since it has tridiagonal form, all the processes with q-matrix A 
will be hirth and death processes. 

By studying A we will be able to identify exactly one process ^t on some 
new probability space (5, 5, V) such that Pn{t) = "Pi^t = n). Since ^ will turn 
out to be dishonest (meaning that "P-a.s. ^ will escape to infinity in finite 
time), the conclusion will be that limj^oo J2T=i'^'^ i-^ni'^)] — ^■ 

5.1. Minimal process. In general, given a g-matrix A, there can be many 
processes x with different laws ynit) = T-'ixit) = ^)i ^H satisfying either 
the forward y' = yA or the backward y' = Ay equations associated with A. 
Whether the solutions of the two systems are unique depends on some well- 
studied properties of the g-matrix. 

In the present case, A is stable (no — oo entries appear in the diagonal) 
and conservative (no mass disappears at zero because //i = 0). It is well 
known that to any stable qi-matrix is associated a process, called minimal, 
whose law satisfies both systems of equations. 

The latter is the naive process that anyone would construct from A, as fol- 
lows. Given a probability space (S, S, V), let ^t be a continuous-time Markov 
chain on the positive integers, with initial distribution 

V{Co = n)=Pn{0), n = l,2,..., 

and jump rates given by A entries, that is, the process waits in a state n 
for an exponential time with rate A„ -|- /i„ and then jumps at n -|- 1 or n — 1 
with probabilities 7r„ and 1 — 7r„, respectively, where 

An 

TTn :~ 



An + A'ti 
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Let r € [0, cxj] denote the first time such that in [0, r) the process has under- 
gone infinitely many jumps. We say that the process reaches the boundary 
at time r and we give no special "return" rule if the process reaches the 
boundary in finite time. Hence, if for uj G S, t{uj) < oo, then Ct{^) is not 
defined for t > t{uj). Notice that, given s > 0, V{t > s) = Yl'^=i'^i^s = n) 
could be less than 1. 

If the minimal solution of a g- matrix is honest, it is the unique solution 
for each one of the two systems and the g-matrix itself is called regular. As 
anticipated, the minimal solution, which is the law of the process described 
above, is not regular if the coefficients kn grow too fast (Proposition 8 below), 
nevertheless it is the unique solution of the forward equations (Proposition 9 
below), while the backward equations have infinite solutions. 

This uniqueness is very important because it ensures that V{£,t = n) = 
Pn{t) := ||x||~^E'3[X^(t)]. If we denote by £ the total energy of X, 



^ oo 

(12) £{t) := -Y^Xlit) 

this means in particular that we can study E*3[£'(t)] through V(t > i), 

^ oo oo 

(13) E«[^(t)] = - J]EQ[x2(i)] = ^(0) 5^p.(t) = ^(O)P(r > t), 

n=l n=l 

which will be the aim of Section 5.2. 

Proposition 8. The q-matrix A is not regular if and only ifYln^^n'^ < °°- 

For the proof we make use of results by Renter and Anderson, which 
are efficiently exposed in the book by the latter [2]. It is, however, not too 
difficult an exercise to prove the "if" direction with elementary notions. 
Truly, Proposition 11 and Lemma 14 below provide such an argument and 
we refer the reader who wants some insight to them. 

Proof of Proposition 8. By Corollary 2.2.5 of [2], in the conservative 
case, the minimal solution is honest if and only if the backward equations 
have a unique solution. 

By Theorem 3.2.2 of [2] the q-matrix of a birth and death process has 

a unique solution of the backward equations if and only if the following 

quantity is infinite: 

oo , 
„ ^ V^ / J_ l^n MnMn-1 /^n • • ' /^2 



^^\^n \i\i-l A„A„_iA„_2 A„---A2Ai 

n=l ^ 

Since An = Mn+i, we get R = Yjn'^K^ = Yjn'^K'^ ■ ^ 
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Proposition 9. The forward system of equations (10), together with 
condition (11) admits a unique solution. 

Here again the proposition can be seen as a simple application of a result 
from the book by Anderson, specifically Theorem 3.2.3 of [2]. 

Uniqueness could also be proved with an analytic approach, based on 
the parabolic structure of the equation which is apparent if we remember 
An = /^n+i and we rewrite (10) as 

(14) -j^Pn{t)=\n{Pn+l{t)-Pn{t))-\n^l{pn{t)-Pn^l{t)), t > 0. 

Nevertheless, we believe it would be interesting to show a completely differ- 
ent and entirely elementary proof that maybe could also be used when the 
parabolic nature is lost and the associated process is no more a simple birth 
and death. 

Proof of Proposition 9. By hnearity we can suppose Pn(0) = 0, with 
condition (11) still holding and —1 <Pn{t) ^ 1- 

Suppose by contradiction that {pn)n&i is a nonzero solution. Without 
loss of generality we can suppose pi(to) = <^ > for some to > [take the 
largest uq such that Pn = for all n < no, so we have Pno{to) > for some 
to > 0; then shift and rename the indexes of the sequence (pn) in such a way 
that the new pi is the old Pno]- Define the partial sums 

n 

(t)n{t):=Y.Pj{t). 
i=i 

We notice that a simple computation starting from (14) yields 

(15) -7:(t>n = K{Pn+l-Pn)- 

Then define the times 

(16) t„:=inf{t|(^„(t) >n(5}, n > 1. 
We claim that for all 7t. > 1, 

tn<tn-l and Pnitn)>S 

so that the sequence {tn)n>o is finite (in fact, decreasing), in contradiction 
with the position (pnit) < 1 for all n and for all t. 

We shall prove the claim by induction. 

For n = 1, by definition ti < to < +oo and pi(ti) = </>i(ti) > 5. 

Let us suppose that the claim holds for n. By the definition of t„, -^(pnitn) > 0. 
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By (15), this implies Pn+i{tn) >Pn{tn) > ^ and hence 

(l>n+l{tn) =(t>nitn) + Pn+litn) > n5 + 6 

thus, tn+i < tn- This imphes that cj)n{tn+i) < n6, so that necessarily 

Pn+l{tn+l) >S. 

The induction and the proof are complete. D 

We remark the fact that given the condition X^^^ ^^n^ < ^^' the forward 
equations have a unique solution while the backward have infinitely many. 
This fact might appear a bit disconcerting if one notices that A is symmet- 
ric and hence, forward and backward equations are formally identical. The 
explanation is that any proper solution Pn{t) of the forward system of equa- 
tions must be summable in the sense that X^„Pn(0 — 1- On the contrary, 
if {qn{t;k)}n is a solution of the backward equations with initial condition 
qn{0', k) = 6k^ny it must satisfy ^^ Qnit', k) <1 for all n. 

5.2. Time of escape. In this section we study the law of r, the time of 
escape to infinity of the minimal process. The main result is Proposition 11, 
which is generalized by Lemma 14. 

Lemma 10. Suppose ^'^i k^ < oo and that the minimal process starts 
from 1 . For n>l, the number of times the minimal process visits state n is 
a geometric r.v. with mean {k'^ + k'^_i)'^'^j^k~ . 

Proof. We follow ideas from Feller [14]. Let pij denote the transition 
probabilities of the discrete time Markov chain embedded in continuous- 
time minimal process and let a^"^' = {an }n>« denote the probabilities that 
the chain starting from states n larger than i will never get to i. Then cr'*^ 
is the maximal solution of 

(17) Xn = ^^Pn,jXj, n>i, 

satisfying < a;„ < 1 for all n. (This solution can be zero.) 

If we let Xi = for sake of notation, in our case the system (17) reduces 
to 

Xn = T ■ X„_i + - ■ Xn+1, n>l + l, 



yielding 



Xn+1 Xn — , \Xn Xn—lj} TL ^ I -\- I, 

An, 
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and then by induction, for n>i, 



2 



^i _ ^ H 



j=k+l ^ 
n.-l 



^n — Xi+l'^i / J "^ri 



'm=i 



By hypothesis the sums are bounded and hence, the maximal solution is 
obtained by choosing Xj+i such that lim^x^ = 1, that is, 

\ -1 



oo 



^» -\ U-i\^ 1,-2 



^i 2_^ " 



'i+l 

\ m=i / 

hence, the chain is transient. 

Now suppose that the chain is starting from 1. It will visit i at least 
once. When it does, the probability that it is the last visit is Pi,i+i(Tl_^i, so 
by strong Markov property, the total number of visits to i is a geometric 
random variable with mean 



oc 
ii) \~^ -(1,2 ^1,2 \Sr 1,-2 

m 



m=i 

as required. D 

Proposition 11. Suppose z^oo := S^i ^^n^ < oo and that the minimal 
process starts from 1. Let Tn be the total time the minimal process spends in 
the state n, 

Tn:=C{t>{):it = n}, 

so that the time of escape at infinity is r = ^^q-^"- 

Then for all n>l, T^ is an exponential r.v. with mean f„ := X]£n ^i" 
and in particular 

oo oo 

n=l n=l 

Moreover, there exists h> such that for all t 

g-iM <V{T>t)<e-^^''^+^. 

Proof. The total time spent in a state n is the sum of many i.i.d. 
exponential waiting times of rate /c^ + k'^^^. 

Since the sum of a geometric number of i.i.d. exponential r.v.'s is expo- 
nential, by Lemma 10, T„ is exponential and its mean is as required. 

The lower bound comes easily from r = X^„T„, since 

V{T>t)>V{Ti>t) = e-^/''\ 
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For the upper bound, let H := — Yl^=i ^n^ogUn and let 

oo 

h := H/voo + log i^oo = - y^ ^^ log 

n=l 

We need to prove that 

lit is such that z^ooc'^"*'' ''°° > 1, we are done. Otherwise, define the sequence 
of numbers {9n)n>i in such a way that for all n, 

The numbers On are positive since Vn ^ z^oo moreover, 

oo oo 



n=l n=l 



-UnlogUn- -{H -t) 

t t ^Cd 



Now, 



(18) 



V{T>t)<v{\J {Tn > Bnt} < J2 'PiTn > Out) 
\n=l / n=l 

oo oo 



n=l n=l 

and the proof is complete. D 

Remark 12. Our argument did not make use of the joint law of the r„'s 
which is unknown. If the r.v.'s r„ were independent, a standard exponential 
bound would yield ^(r > t) < e~*/'^'^+^ ^ so we have the strong feeling that 
1/z^oo is not a sharp bound for the true rate, which could actually be l/i^i. 

It should also be noted that one cannot get rid of h in the previous 
statement; one can prove that ^log7'(r > i)|i=o = 0. 

The following lemma clarifies that either r = oo a.s. or r belongs to any 
interval [a, 6] with positive probability. 

Lemma 13. Suppose V(t > T) < 1 for some T. Then V{t > a) > V{t > b) 
whenever a <b and in particular V{t > t) < 1 for all t> 0. 

Proof. For i > 1 and t > 0, let 

oo 

qi{t) := V{t > t\Co =i)= Y,Pi,n{t)- 

n=l 

By Chapman-Kolmogorov, 

oo 
(19) V{t >t)=Y, Pmit - S)qm{s) <V{T>t- s). 

m=l 
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Suppose by contradiction that Vij > a) = V{t > b), then letting s = b — a, 
t = b,we have equahty in (19). This imphes that qm{s) = 1 for all m so that 
we get equality for any choice of t > 0, meaning that V{t > •) is periodic as 
well as nonincreasing. The only possibility is V{t > t) = 1 for all t. D 

We will need the following lemma, that is, the probability of no explosion 
is larger if one starts from 1 than in any other case. This is a well-known 
fact whose proof is a standard exercise we do not repeat here. 

Lemma 14. V{t >t)< V{t > t|^o = 1) ■ 

6. Bounds on the energy. We are finally able to give statements on the 
decay of the energy 8{t) as i — > oo. First of all we obtain exponential esti- 
mates under Q, both in C^ and pathwise. Then we introduce a smallness 
condition on £{0) which is what we need to translate these results under P. 
The section concludes with the proof of Theorem 1 . 

Proposition 15. Suppose Voo '■= S^i "-^n^ < ^o. Let X be an energy 
controlled solution and denote by £{t) := 2X]n^i^n(0 ^^-^ total energy at 
time t. Then 

limEQ[£'(t)]=0. 

In particular there exists a positive number h such that for all t>0 
(20) EQ[^(t)]<e-*/'^-+^.S(0). 

Proof. By definition 

E^[^(i)] = lj2^'^^^nit)]=£{0)Y.Pn{t), 
n n 

whence, by Proposition 9, 

= £{0) Y, nit = n) = £{Q>)V{it < cx)) = f (O)P(r > t). 

n 

Finally, we apply Lemma 14 and Proposition 11 and we find (20). D 

Using Borel-Cantelli arguments, one can deduce from (20) some Q-a.s. 
statements about the decrease to zero of the energy, at least on given se- 
quences of times going to infinity. To extend this to all sequences we will 
need the following lemma. 

Lemma 16. Let X be an energy controlled solution and denote by £{t) := 
2 Sn^i ^n(^) ^^-^ total energy at time t. Then for every t> s>0 we have 

Q{£{t)<£{s)) = l. 
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Proof. Let s > be given and set x = -^{^)- Consider the linear equa- 
tion on [s, oo) with initial condition x- It is proved in [4] that it has a unique 
strong solution y, with the property 

(oo oo \ 

n=l n=l / 

for every t G [s, c») (the result in [4] is for constant initial conditions, but 
the extension to nonanticipative square integrable random initial conditions 
is straightforward). But also X restricted to [s,cx)) is a solution of the same 
equation, hence, equal to Y, on [s,cxd). The previous identity is thus equal 
to the claim of the lemma. D 

Now we can prove Q-a.s. exponential decay of energy. 

Proposition 17. Under the same hypothesis of Proposition 15, the total 
energy of solutions goes to zero at least exponentially fast pathwise under Q, 

limsup — logiS(i) < , Q-a.s. 

i— >oo t Uqq 

Proof. Let e > be given. Set a := I/z/qo + e. We have 

Qin~Hog£{n) > a) < e-""E«[f (n)] < Ce"^'^, 

where, by Proposition 15, C = e £{0) does not depend on n. Hence, the 
above probabilities are summable on n and by Borel-Cantelli lemma there 
exists a measurable set A'^ with Q{N) = and the following property: for 
every uj G N'^ there exists nQ{uj) such that, for all n > uq^uj), £{n,uj) < e""". 
Taking the supremum for n < no(c<j), we obtain that there exists a constant 
C{uj) > such that 

£{n,uj)<C{uj)e-"'' 

for all a; G N^ and n > 0. From Lemma 16, there exists a^measurable set N 
with Q{N) = such that £{r,uj) < £{\r\,uj) for all lo G N'^ and all rational 
numbers r G [0, oo). This implies, for uj G A^'^ Pi A^*^, 

£{r,uj) < £{[r\,uj) < C(6j)e-"W < C'(w)e~"" 

with C'{uj) = C7(w)e"', for ah r G [0,oo) nQ. 

With Q probability one, the function £{t) is lower semicontinuous, being 
the supremum in A'^ of the functions X^„=i ^n(^) which are continuous. Thus 
we get 

Q{£{t) < C'e-"* for every t > 0) = 1. 

Letting e go to zero on the rationals completes the proof. D 
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The reader should be aware that this proposition does not automatically 
hold under P. Truly, P and Q are equivalent on all Ft, but C" is Fqo- 
measurable and not i<t-nieasurable for any t. 

In general, we cannot prove that P and Q are equivalent on Foq, so we 
cannot translate such claim into a similar statement on the original nonlinear 
equation (1). 

When £{0) is small enough, however, we can prove the equivalence and 
hence, compute exponential upper bounds for S{t), both P-a.s. and in mean 
value. 

Proposition 18. Let X be an energy controlled solution and denote by 
£{t) := i X^^x -^ni'^) ^^-^ total energy at time t. Suppose v^o '■= S^i "-^n ^ < 
oo and z^oo^(O) < 1. Then 

so that in particular, P and Q are equivalent on Fao- 

Proof. By Proposition 15 and the definition of energy controlled solu- 
tion, we have 

E'^[£{t)] < ^(0)6"*/''°°+'' Vt > 0, 

0<£{t)<£{0), Q-a.s. 

Let X := f^ £{t) dt > and x > 0. Then 

xQ{X>x)<E^[X-X>x\= / E^[£{t);X>x\dt 

Jo 

POO 

< / mm(E^[£{t)]-£{0)Q{X>x))dt 
Jo 

POD 

<£{0) / min(e-*/''-+^;Q(X>x))dt 
Jo 

= £{0) Q{X>x)dt + £{0) e-^'^^^^dt, 

Jo Ju 

where u is such that e~"/''°°+'' = Q{X > x). Hence, 

xQ{X >x)< £{0)Q{X > x)u + 5(0) 7^006""/''°°+'' = £{0)Q{X >x){u + v^c)- 

If Q{X > x) = for some x > 0, then clearly X is bounded and we are done. 
Otherwise we get 

that is 

Q{X > x) = 6""/''°°+'' < e.-^/i'^^^m+h+i^ 
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yielding 

and finally 

/"CO /"OO 

K^[e^] = / Q{e-'^ >y)dy<l + e^+^ / y-i/(^oo£(o)) ^^ ^ ^^ 
where we used i^oo^(O) < 1. □ 

The heuristic behind the proof, which is quite hidden, that is, given that 
E'^[X] is bounded, E'^[e ] is maximum if X is spread as much as possible. 
Since X = J £{t)dt and £{t) £ [0,£{0)], this is done by choosing £(t,uj) € 
{0,£'(0)}, and in particular £{t,uj) = £'(0)I[o,j/(a;)](t). 

Corollary 19. Under the same hypothesis of Proposition 18, meaning 
in particular that £"(0) < 1/v^, the energy goes to zero at least exponentially 
fast pathwise under P, 

limsup-log<S(t) < , P-a.s. 

Proof. Just a direct consequence of Propositions 17 and 18. D 

The same condition on the smallness of £'(0) arises when we want to 
establish an exponential decay for the mean value of £{t) under P. 

Proposition 20. // 

then 

(21) E^[£{t)]<£{0)exp{{l-l/p)[h + {p£{0)-a)t]) 

for every p> 1. In particular, under the same hypothesis of Proposition 18, 



1, „p . 1 



2 



(22) limsup-logE^[6'(t)] < (1 - \/£{'d)uoo) 

t->oo t l^oo 

Proof. The density ft of P with respect to Q on Ft is 

(23) /i = exp(Mj + i[M]i), 
where 

OO „t Pt °° 

Mt:=y2 Xnis)dWn{s), [M]t= / Vx2(s)ds. 

n=l'^0 Jo „=i 
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From (6) we have 

(24) exp(A[M]i) < e"^^^^^^' 

for every A > 0. 

For every p,p' > 1 with - + — = 1, from the a.s. condition £(t) < £{0) and 
the assumption of the proposition we have 

E^[£{t)] =E'^[ft£{t)] <KQ[fj']^/PEQ[£{tf']'^/p' 

< E^[f[]'^^PE^[£{t)£{Of'~^]^/P' 

From (23) we have 



:E-' 



:E-P 



exp{{p-l)Mt + ^^-^[M]t 



exp( ip-l)Mt - (^^[A/]i] exp(^^^[M]t 



and now we use (24) to get 

< g(p-l)p^(0)i]gP[g(p-l)Mt-((p-l)V2)[M]tj ^ g(p-l)p£(0)t ^ g(pVp')£(0)t 

by Girsanov's theorem. To summarize: 

EP[£{t)] < £{0)e^yp')pmte-Hp')t+h/p' 

which imphes the first claim of the proposition. 

To prove the last statement, let a = l/i'oo- Then optimization on p under 
the condition £{0)l'oo < 1 gives that the right-hand side of (21) is mini- 
mum when p is equal to (f>{t) = \/ ftn) ■ [We notice that (/)(t) > 1 for 



t > , / „<7/n\ ■] With a simple computation, (21) becomes 



E^[£{t)]<£{0)expi-{p-l)(^-^-l]£{0)t 
Letting p = <^(t) > 1, we get 

ilogEnSWI < ilogf (0) - (m - Ifm). '> yj'_n,y 

Taking the limsup for t —^ oo leads to (22). D 
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6.1. Proof of Theorem 1. We have u^ = X^^i^^n^ < oo, so by virtue 
of Proposition 8, V{t > i) < 1 for some t. Then by Lemma 13, V{t > i) < 1 
for all t. By equation (13) ¥P[£{t)] < £{0) and since £{t) < £{0) Q-a.s. we 
get Q{£(t) = £{0)) < 1, for all t. Equivalence of P and Q on Ft yields the 
first statement. 

By Proposition 15, Q{£{t) > e) < Ke~^''^°° for some K not depending 
on t, hence, for t large enough, this event has Q- (and hence P-) probability 
less than 1, so we proved the second statement. 

The third statement is proved in Corollary 19 and Proposition 20. 

The proof of Theorem 1 is complete. 

7. Lack of regular solutions. As a consequence of our result on the dis- 
sipation of energy we can prove that there exists no regular solution. 
Define the space 

which is an Hilbert space under the norm ||a;||y = X^^x ^n^n- 

Proposition 21. Assume that {X{t);t G [0,r]} is an energy controlled 
solution. Then, 

p( I ||X(t)||^dt = oo] >0. 

Proof. We will actually prove the following statement. Assume that 
{X{t);t G [0,r]} is an energy controlled solution such that 

(25) p( I ||X(t)||^dt<ooj =1. 

Then, for every t G [0,r], P{£it) = £{0)) = 1. 

Since the latter is in contradiction with Theorem 1, then (25) will be 
proven to be false. 

Let X be a solution as in the claim. By Ito's formula, 

(N \ N 

Y,Xl =2^{kn^iXl_,Xn - KXlXn+l) dt 
n=l / n=l 

N N 

- Y.^kl + kl_^)Xl dt + ^(fctiX2„i + klxl^^) dt 

71=1 n=l 

N 
+ 2 ^{kn-lXn-lXn dWn-1 - A;„X„X„+i dWn) 



n=l 
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^ —2ki\iXj^Xj\[^i dt 

— 2kNX]\[X]\!^i dWN- 
Hence, 

TV N t 

^Xlit)-^xl= {-2kNXlXN+i + klXl^, - k%Xl) ds 

n=l n=l -^0 

(26) ^ 

2kNXNXN+ldWN{s). 



Jo 



'0 

Notice that |X7v+i| < ||3;||;2, k]\f > 1 and kjy < k^j^i, so that the first integral 
can be bounded by 

C I {kj^X^ + k^^-^^Xjy^i + k^Xj^) ds. 
Jo 

Then the a.s. inequahty j^ \\X{t)\\ydt < cx) imphes that Jq k'j^X'j^ds goes 
to zero a.s. and in probabihty. 

The latter is true also for the stochastic integral in (26) by continuity in 
probability of stochastic integrals; since k^X^X^j^i converges to zero in 
probability in L^(0, T), its stochastic integral converges to zero in probabil- 
ity. This can be easily proved by applying the following well-known inequality 
of stochastic integrals. For every 7, 5 > 0, 



P 



t 

kNXNXN+ldWN{s) 




>7)<P( / k%X'^X'^^,ds>s]+^. 



I klxlxl^^ ds>5\ + 



Finally, all terms on the RHS of (26) converge to zero in probability. We 
get 

N 

,.9..s P ,, , 

\x\ 



p. Il„l|2 
n=l 



N-^od" "^^' 



n=l 

But we know, by P-a.s. monotonicity, that 

N 00 

E^n(t) AE^n(t). 
n=l n=l 

Hence, the claim is proved and the proposition as well. D 



Remark 22. If one could prove local existence (up to some random 
time) of solutions in V with initial conditions in V, the above proposition 
would show a blow-up in the V norm. This appears to be a difficult open 
problem in the stochastic case, while in the deterministic case, it is well 
known (see, e.g., [7]). 
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